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^ i Abstract 

Ph ■ An eikonalized elastic proton-proton and proton-antiproton scat- 

Q\\ tering amplitude F(s,t), based on the suggestion of a finite sum of 

ladder diagrams, calculated from QCD and the number of s— channel 
gluon rungs and correspondingly the powers of logarithms in total 
| cross section depends on available increasing energy. Explicit expres- 

sions for total cross section involving three and four rungs (four and 
^4 five prongs) with In 3 (s) and Zn 4 (s) as highest terms, respectively) are 

r~| ! fitted to the all available proton-proton and proton-antiproton total 

cross section data. Predictions for pp total cross section at LHC en- 
ergy are given and compared with prediction of several Regge-models 
including possible hard Pomeron contribution. 

1— I ■ 

1 Introduction 

The renewed interest in high-energy soft interactions is largely triggered by 
the recent LHC measurements by the TOTEM pQ measurements of pp scat- 
tering. The accessible energies offer new possibilities to check the basic theo- 
retical concepts, such as unitarity and the properties of the (multy) Pomeron 
exchanges. 

In the present paper we develop further and calculate the effects coming 
from the inclusion of both a "soft" and "hard" Pomeron, the Pomeron being 
a finite sum of reggeized gluon ladders, different from the " QCD Pomeron" , 
which is an infinite sum of such ladders [21 El H] , resulting in the so-called 
supercritical behavior of the total cross section. In that approach, the main 
contribution to the inelastic amplitude and to the absorptive part of the elas- 
tic amplitude in the forward direction arises from the multi-Regge kinematics 
in the limit s — > oo in leading logarithmic approximation. In the next-to- 
leading logarithmic approximation, corrections require also the contribution 
from the quasi- multi-Regge kinematics [3]. Hence, the subenergies between 
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neighboring s-channel gluons must be large enough to remain in the Regge 
domain. At finite total energies, this implies that the amplitude is repre- 
sented by a finite sum of iV terms [6J, where N increases like Ins, rather 
than by the solution of the BFKL integral equation [2j [3j 0]. The interest 
in the first few terms of the series is related to the fact that the energies 
reached by the present accelerators are not high enough to accommodate for 
a large number of s-channel gluons that eventually hadronize and give rise 
to clusters of secondary particles [7]. Consequently, one can expand the 
"supercritical" Pomeron ~ in powers of ln(s). 

In Ref. [6J a model for the Pomeron at t = based on the idea of a finite 
sum of ladder diagrams in QCD was suggested. According to the idea of that 
paper, the number of s-channel gluon rungs and correspondingly the powers 
of logarithms in the forward scattering amplitude depends on the phase space 
(energy) available, i.e. as energy increases, progressively new prongs with 
additional gluon rungs in the s-channel open. Explicit expressions for the 
total cross section involving two and three rungs or, alternatively, three and 
four prongs (with ln 2 (s) and ln 3 (s) as highest terms, respectively) were fitted 
to the proton-proton and proton-antiproton total cross section data in the 
accelerator region. In a related paper [8] the Pomeron was considered as a 
finite series of ladder diagrams, including one gluon rung besides the Low- 
Nussinov " Born term" [9] and resulting in a constant plus logarithmic term 
in the total cross section. With a sub-leading Regge term added, good fits 
to pp and pp total as well as differential cross section were obtained in [8]. 
There is considered the opening channels (in s) as threshold effects, the 
relevant prongs being separated in rapidity by In Sq, sq being a parameter 
related to the average subenergy in the ladder. Within the "finite gluon 
ladder approach" to the Pomeron (see [BJ and references therein), several 
options are possible. In Ref. [5] a system of interconnected equations was 
solved with several free parameters, including the value of Sq, that determine 
the opening of each threshold (prong). In that paper finite gluon ladders 
were calculated from QCD, where the important dynamical information is 
contained in p of Eq. (7) of that paper, including In s terms multiplied by the 
QCD running constant a s constraining the interconnection between various 
powers of the logarithms in the total cross section. In [10] an unitarization 
procedure was also included: the QCD-inspired amplitude was treated as a 
Born term, subject to a subsequent unitarization procedure. 

In present paper we use the abovementioned "finite gluon ladder ap- 
proach" to the unitarized QCD-inspired amplitude of [TO] and calculate the 
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values of pp total cross section at LHC energy 7 and 14 TeV for several 
possible So- For the sake of completeness, we compare the pp total cross 
section resulting from this model with those from popular alternatives, such 
as: exponential, double logarithmic form and the DL two-Pomeron model 
[TT] . Fits to the pp and pp data is performed from yfs = 5GeV up to highest 
energy data PQ, [T2] . 

2 Total cross sections from a finite sum of 
gluon ladders 

Following [SI, we write the Pomeron contribution to the total cross section 
in the form 

JV 

a ( t p \s)^j:f i e(s-si)e(s^-s), (1) 

t=0 

where 

/< = £ a l3 U , (2) 

3=0 

So is the prong threshold, 9{x) is the step function and L = ln(s). Here and 
in the following, by s and s respectively, s/(lGeV 2 ) and s /(lGeV 2 ) are 
implied. The main assumption in Eq. (QQ) is that the widths of the rapidity 
gaps ln(so) are the same along the ladder. The functions fi are polynomials 
in L of degree i, corresponding to finite gluon ladder diagrams in QCD, 
where each power of the logarithm collects all the relevant diagrams. When 
s increases and reaches a new threshold, a new prong opens adding a new 
power in L. In the energy region between two neighboring thresholds, the 
corresponding fi, given in Eq. (pp), is supposed to represent adequately the 
total cross section. 

In Eq. ([1]) the sum over N is a finite one, since N is proportional to 
ln(s), where s is the present squared cm. energy. Hence, this model is quite 
different from the "canonical" approach [2], where, in the limit s — > oo, the 
infinite sum of the leading logarithmic contributions gives rise to an integral 
equation for the amplitude. 

To make the idea clear, we first describe the mechanism in the case of 
three gaps (two rungs). In this case 

fo(s) = a 00 (3) 
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for 



s < s , 
fi(s) = a w + a n L 



(4) 
(5) 



for 



So<S< Sq, 



(6) 



/ 2 (s) = a 2 o + «2i-^ + a 2 2-^ 2 



(7) 
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By imposing the requirement of continuity (of the cross section and of its 
first derivative) one constrains the parameters. 

The same procedure can be repeated for any number of gaps. 

Notice that the values of the parameters depend on the energy range of 
the fitting procedure. For example, the values of the parameters in f is 
fitted in "its" range, i.e. for s < sq, will get modified in fi with the higher 
energy data and correspondingly higher order diagrams included. 

As a first attempt, only three rapidity gaps, that correspond to two gluon 
rungs in the ladder were considered [5J. Fits to the pp and pp data were 
performed up to the highest energy Tevatron data 1.8 TeV. The value of the 
rapidity gap turned out to be ^/so ~ 12 GeV, i.e. the value for which the 
energy range considered is covered with equal rapidity gaps uniformly. 

3 Fit with the Unitarized QCD-inspired Pomeron 

From the theoretical point of view, the phenomenological model of Section 
2 corresponds to the explicit evaluation in QCD of gluonic ladders with an 
increasing number of s-channel gluons. This correspondence is far from literal 
since each term of the BFKL series takes into account only the dominant 
logarithm in the limit s — > oo. In the following we give concrete expressions 
for the forward high energy scattering amplitudes for hadrons in the form 
of an expansion in powers of large logarithms in the leading logarithmic 
approximation. 

We start from known results obtained in paper [3] where an explicit ex- 
pression for the total cross section for hadron-hadron scattering has been 
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obtained. In the high energy limit, it is convenient to introduce the Mellin 
transform of the amplitude 
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To obtain the total cross section of proton-proton scattering it was used 
[TO] the ansatz of Ref.[13] for the impact factor of a hadron in terms of its 
form factor $(g 2 ): 

$o(fc) = ak 2 e~ bk2 , (10) 

as a result 
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where p defined as 



^(ln2)+4(ln2) 3 -^C(3) 



3a, 



7T 



In s 



(11) 
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and C(3) is the Riemann's Zeta function £(3) ~ 1.202. 

This approach was used [10] to predict the total cross section and ratio of 
real part of elastic amplitude to it imaginary part at LHC energy. The strong 
coupling a s is assumed to be frozen at a suitable scale set, for example, by 
the external particles. 

Now we proceed with a fit and predictions with another Pomeron contri- 
bution corresponding to four rungs: 



a t = (7ra7(26))[l + Cl p + c 2 p 2 + (c 3 /3)p 3 + (c 4 /12)p 4 + ...], 



(13) 



where 



ci = 2(ln2), 



c 2 = ^ + 2(ln2) 



c 3 = 3,8511 and c 4 = 10,145. 
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Apart the calculations of the total cross section involving fits of a vari- 
ety of Pomeron models, it is interesting to estimate the value of Sq, a basic 
parameter in the finite series of QCD diagrams. Each set is 'active' in its 
rapidity gap, i.e. the parameters in (jSJ) should be fitted in each energy 
interval separately and relevant solutions should match by imposing conti- 
nuity of the total cross section and of its first derivative. Here we present 
the result of a fits to the existing experimental data (including the new LHC 
data [1]) and predict the values of the total cross section at the next LHC 
energy with the unitarized contribution of the explicit BFKL iteration. As 
will be seen below, the values of so are quite constrained by the fits. 

First of all, in the framework of our approach we mast put ^/so > 5 GeV. 
As a convenient minimal value of rapidity gap we chosen ^/sq = 6.8 GeV, 
where in estimating <Jtot{pp) we are restricted by the maximal power Zn 4 (s). 
From the experimental data y/so < 24 GeV, one expects a high-quality 
description of the total cross section with the Pomeron contribution of Eq. 
(3), including relevant secondary Regge-pole terms. To cover the whole fitted 
range, we choose the interval 10.9 GeV < ^/so < 24 GeV, where four gaps 
are sufficient to estimate of total cross section at the maximal LHC energy. 

h{s) = a 30 + a 31 L + a 32 L 2 + a 33 L 3 (14) 

for 

s 3 < s < 4 (15) 

For the remaining part of the studied region 6.8 GeV < ^/sq < 10.9 GeV it 
is necessary to take into account five gaps to estimate the total cross section 
at the LHC. 

fi(s) = a 40 + a 41 L + a i2 L 2 + a 43 L 3 + a 44 L 4 . (16) 

for 

s 4 <s< 4 (17) 

Below we perform such calculations in the framework of the eikonal for- 
malism and compare the results with the experimental data. We start from 
Eq. (fHfl) for the pp and pp total cross section. Supplying that expressing with 
an exponential t-dependence, we get the elastic scattering amplitude: 

F Born (s,t) = A(-i~s) 1+a ' t x 

[oq + ai7 ln(-is) + a 2 7 2 \n 2 (-is) + a 3 >y 3 ln 3 (-is)}e Bt , (18) 
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where a' and B are fitting parameters, 7 determined as (j!2j) and 
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y ln2 + 41n d 2- i C(3) 



A= — 
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(19) 



In the eikonalization procedure we follow Ref. [H], according to which 
the Pomeron amplitude 



F P (s,t) = is J™ 6d&J (V=t)(l - e ix{b ' s) 
where Jo is the Bessel function of zeroth order and the eikonal x is 



s Jo 



(20) 



(21) 



Inserting the expression for the Pomeron into Eq. f[2~Tj) and expanding 
the exponential in (l2"0l . one find for the eikonalized Pomeron amplitude 



The corresponding forward Pomeron amplitude is 

F P (s, t = 0)= 2isfi[C + ln{£/n) + £1 (£///)], 

where 

fj, = B + a'ln(— is); 
A 



and 



£0 = o Q , & = ai7ln(-is) 



(22) 

(23) 

(24) 
(25) 
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Figure 1 : Total pp ( upper curve ) and pp cross sections from the uniratized ( eikon- 
alized) version of the model. The vertical lines corresponds to boundaries of rapidity 
gaps. 



6 = a 27 2 ln 2 HS), £3 = a 3 7 3 ln 3 (-*S), (26) 

where C = 0.577216 is the Euler constant and E\ is the asymptotic form of 
the first order exponential integral: 
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For the case of four rungs, the Born amplitude is 



+ 



where 



F B om{s,t) = A{-i§) 1+a>t [a + ai7ln(-is) + a 2 ^ 2 \n 2 {-iS) + 
+ a 3 >y 3 ln 3 (-is) + a 4 <y 4 ln 4 (-is)]e B \ 



(27) 



(28) 



7T 



2 -7T 2 



} ln2 + 41n 3 2- -C(3)j7 2 + 2/n2. 



37T 2 2 7T 2 nJ „ 

a 2 = — c 4 7 + yjj + 2/n 2 2, 



a 3 



lrvr 2 



ln2 + 41n J 2- -C(3) 
2 4 



° 4 = 12 ° 47 



where A, \i and 7 are the same. 

A 



and 



£= 2^2(£o + 6 + 6 + 6 + &), 



£0 = o ; 6 = aaln(-is); £ 2 = «27 In (-is); 



(29) 



£3 = a 3 7 3 ln 3 (-is); £4 = a 4 7 4 hi 4 (-is). 

The obtained eikonalized Pomeron terms are appended by a contributions 
from secondary reggeons, p and ui: 

F±(s,t = 0) = g f s a f (0) ±i goj s a " io \ (30) 

where the +(— ) sign corresponds to pp(pp) scattering, the resulting forward 
amplitude being 

F%(s,t = 0) = F P (s,t = 0) + F£(s,t = 0). (31) 
For the total cross section the norm 

a = —ImFPP(s,t = 0) (32) 

was used. 

Consider in more detail the fitting procedure in the approach with dif- 
ferent numbers of rungs (power of log's). In the case of 3 rungs, for values 
^/Iq within the interval 10.9 GeV - 24 GeV the whole range of the data 
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5 GeV < y/s is covered. We have chosen the parameters B = 0.116 and 
a' = 0.134 [TO]. 

The parameters of secondary reggeons are the same for whole fitted ex- 
perimental data region. Notice that in this approach for the Pomeron contri- 
bution there is only one free parameter in every separate interval. Therefore, 
for its determination, in principle, it is sufficient to have one experimental 
point. In our case this condition is fulfilled. 

To estimate So in the remaining part of the whole investigated region 
6, 8 GeV < t/so < 10.9 GeV, we added the new LHC [lj data as well as 
the less reliable cosmic-ray data [15]. In our calculations so is not a free 
parameter. Instead, we have performed a series of fits and have chosen a 
value s for the cases where x 2 /dof ~ 1. We find, that this condition is valid 
in the whole interval 6.8 GeV < y/so < 10.9 GeV , where a s is rather small 
~ 0.08 and the prediction of total cross section at 7 TeV LHC energy crosses 
the experimental value within the error bars. In case three rungs, however, 
the predicted value is unacceptable low. In table 1. and fig 1. two of the 
best fits for s = 49 GeV 2 and s = 121 GeV 2 are quoted. Throughout this 
paper we chosen the fits with (x 2 /dof) < 1.0. 

4 (7 to t{s) fits for different Pomeron models 

Here we repeat our previous fit pi)] limiting ourselves by 163 data points 
on pp and pp total cross-sections in the range 5 GeV < \/s < 1.8 TeV 
[12], however adding the new LHC experimental point at 7 TeV with the 
results and predictions for pp total cross section at the LHC energies quoted 
in Table 2. For the completeness and comparison we performed also the fits 
of known Regge-models with different the Pomerons used. The variety of the 
Pomeron- mo dels [1_6] for the total cross section, can be distributed to three 
groups: 

Tripole Pomeron 

ImF p (s, t = 0) = ao + ailn 2 (s) , (33) 

Supercritical soft Pomeron 

ImFp (s, t = 0) = a + ai S aa . (34) 

To this end we have added the two-component Donnachie-Landshoff Pomeron 
model [IT] 

ImFp (s, t = 0) = a + + a 2 s ah . (35) 
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Table 1: Values of the fitted parameters and calculated total cross sections 
at the LHC energies. A in the (18) for different % rung 





3 rungs 




4 rungs 




Parameter 


value 


error 


value 


error 


A, 


3.61 


1.11 


4.83 


0.25 


A l 


1.77 


0.92 


3.44 


0.25 


A 2 


1.06 


0.64 


2.63 


0.26 


A, 


0.59 


0.46 


2.09 


0.27 


A A 






1.66 


0.27 


a s 


0.145 


0.049 


0.084 


0.007 


9f 


-6.17 


0.53 


-7.67 


0.32 


«/(0) 


0.681 


0.028 


0.633 


0.011 


9u 


3.61 


0.17 


3.72 


0.18 


a w (0) 


0.463 


0.015 


0.452 


0.015 


s , GeV 2 


121.0 




49.0 




a(pp),mb 










7TeV 


92.2 


1.5 


94.8 


1.5 


14 TeV 


102.6 


2.1 


106.6 


1.9 
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The hard Pomeron contribution to the pp- and pp- total cross section 
is much smaller than the soft one and it is hard to determine it directly. 
However, as it was shown in [17] ; the data for the total cross section and 
for real part of elastic amplitude indicate the presence of a hard Pomeron 
in up and Kp elastic scattering at t — 0, compatible with that observed in 
deep inelastic scattering. The inclusion of these two Pomerons, together with 
use of the integral dispersion relation and subleading Reggeons, leads to a 
successful description of all hadron and photon scattering on proton data for 
5 GeV < s < 100 GeV [17]. The two- Pomeron approach is also compatible 
with pp and pp data, provided one unitarizes it at high energies. It was found 
that the contribution of the hard Pomeron relative to the soft one is 0.17/55 
(see Tab. 4. in [17j). In the last case we eikonalized the Pomeron terms (see 
the previous part). 

For the soft and hard Pomeron, the eikonalized contribution has the form 

£ = & + & + 6, 

with 

£o = flo, 6 = ai (~s) 5s , 
corresponding to to the soft Pomeron and 

6 = a 2 (-s) S » 

corresponding to the hard one, where 5 S = 0.08 and Sh = 0.45 [17]. To 
perform the fit we need to include the contributions from secondary reggeons, 
p and uj (150"j) . The normalization of the total cross section is the same as in 
the previous secttion (1521) . Fits of all models to pp and pp total cross section 
data were performed in the range 5 GeV < y/s < 7 TeV, the quality of fit 
being x 2 /dof < 1.0. The results along with the expected values for the total 
cross section at maximal LHC energy are quoted in Table 2. 

All calculated values fit the new experimental data p] and predict almost 
the same values at next LHC energy except the unitarized two-pomeron 
model, which predict the value exceeding 120 mb. 
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Table 2: Prediction of Pomeron models for total cross section at LHC energies 





&tot (pp) , mb 






7 


14 


tripole 


98.4 ± 0.9 


110.4± 0.7 


soft 


98.7 ±0.9 


111.8±1.0 


sof+hard 


101.8 ± 1.4 


120.1 ± 1.7 


QCD, 2 rungs 


96.5 ±0.5 


107.6 ± 0.6 


QCD, 3 rungs 


94.9 ±1.3 


105.7 ± 1.3 


QCD, 4 rungs 


96.8 ± 1.3 


107.1 ± 1.5 



5 Conclusions 

Our main goal was an adequate picture of the Pomeron exchange at t = 0. 
In our opinion, it is neither an infinite sum of gluon ladders as in the BFKL 
approach [2J [3j H], nor its power expansion. In fact, the finite series - call 
it "threshold approach" - considered in Sec. 2. and 4. and in the previous 
papers [6], realizes a non-trivial dynamical balance between the total reaction 
energy and the subenergies equally partitioned between the multiperipheral 
ladders. 

The role and the value of the width of the gap Sq, is an important physical 
parameter, independent of the model presented above. We have fitted it and 
obtained the value which exceeds our previous estimates [6] . Another goal of 
the present investigation was the prediction of the proton-proton total cross 
section at the maximal LHC energy with several Regge- and QCD-inspired 
Pomeron models. We checked the traditional Regge-models as well as the 
controversial case of the hard Pomeron. In this case we predict the largest 
value for the total cross section, f» 120 mb (see also [IE]). All model fits 
give values of a (s) compatible with the experiment at 7 TeV within the error 
bars and consequently cannot be ruled out. 
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